Photoinduced melting of charge order in a quarter-filled electron system coupled with 
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Photoinduced melting of charge order is calculated by using the exact many-electron wave function 
coupled with classically treated phonons in the one-dimensional quarter-filled Hubbard model with 
Peierls and Holstein types of electron-phonon couplings. The model parameters are taken from recent 
experiments on (EDO-TTF)2PF6 (EDO-TTF=ethylenedioxy-tetrathiafulvalene) with (QUO) charge 
order, where transfer integrals are modulated by molecular displacements (bond-coupled phonons) 
and site energies by molecular deformations (charge-coupled phonons). The charge-transfer pho- 
toexcitation from (0110) to (0200) configurations and that from (0110) to (1010) configurations have 
different energies. The corresponding excited states have different shapes of adiabatic potentials as 
a function of these two phonon amplitudes. The adiabatic potentials are shown to be useful in 
understanding differences in the photoinduced charge dynamics and the efficiency of melting, which 
depend not only on the excitation energy but also on the relative phonon frequency of the bond- 
and charge-coupled phonons. 
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I. INTRODUCTION 

Introduction of electrons and holes by photoirradiation 
of an insulator may trigger a phase transition, called a 
photoinduced phase transition, which are often charac- 
terized by a nonlinear relation between the field and the 
response. The electronic and structural states evolve in 
a cooperative manner, leading to macroscopic changes in 
the physical properties»ii^ They may create novel tran- 
sient electronic states or dynamical functions. Thus, it 
is important to search for possibility of controlling pho- 
toinduced phase transition dynamics. 

When an ordered state is stabilized by different 
types of interactions, e.g., by electron-electron (e-e) and 
electron-phonon (e-ph) interactions or by different types 
of e-ph interactions, their relative contributions would be 
clarified by investigating the photoinduced dynamics of 
the electron-phonon state. For instance, when an e-ph 
interaction stabilizes the order, the motion of photoin- 
duced charge transfers and that of phonons are strongly 
coupled and often show a coherent oscillation in transient 
reflectivity,^i^i^ 

In most cases, one type of phonons has been consid- 
ered in discussing the stabilization of an ordered state 
and its photoinduced dynamics.^''* If different types of 
phonons contribute to forming a charge order by modu- 
lating different electronic parameters, complex dynamics 
are expected and their analysis would give valuable in- 
formation that cannot be easily obtained by measuring 
the physical properties in thermal equilibrium. 

In this paper, we consider the organic salt (EDO- 
TTF)2PF6 with 3/4-filled band (1/4-filled in terms of 
holes), which shows a metal-insulator transition with 
(0110) charge order accompanied by a large structural 
change based on a doubling of unit cells and deformation 
of EDO-TTF molecules at low temperatures,^ The (0110) 



charge order is thus stabilized by two different types 
of phonons, transfer-integral-modulating displacements 
(bond-coupled phonons) and site-energy-modulating de- 
formations (charge-coupled phonons). An ultrafast tran- 
sition (within 3 ps) from the charge-ordered insulator 
phase to the "metal" phase is induced by weak laser light 
with a pulse width of 0.12 ps and of excitation photon 
energy 1.55 eV— nearly resonant with charge transfer to- 
ward the (0200) configuratioUfi^ The observed spectro- 
scopic change shows cooperative melting of the charge 
order assisted by coherent phonon generationi^iii 

For quarter-filled systems, different charge-transfer ex- 
cited states generally have different shapes of adiabatic 
potentials and thus different directions of initial force ap- 
plied to phonon amplitudes. The photoinduced dynamics 
would have respective efficiency, and some of them can 
be sensitive to the relative phonon frequency, as shown 
in this paper. This fact may be useful in identifying the 
phonon responsible for the coherent oscillation, and more 
generally, in designing a material with more efficient tran- 
sition dynamics. 



II. PEIERLS-HOLSTEIN-HUBBARD MODEL 

As a model for (EDO-TTF)2PF6, we treat the one- 
dimensional quarter-filled Hubbard model with Peierls 
and Holstein types of e-ph couplings, which modulate 
transfer integrals and site energies, respectively. 
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where Cj ^ {cj,^) is the creation (annihilation) opera- 
tor of a hole with spin a at site j, rij^a = Cj ^cj^a, 
"^j ~ ^j,T~l""'j,ii % is displacement of the jth molecule 
from equilibrium, Vj is its deformation amplitude, and iij 
and iij are the time derivatives of Uj and Vj , respectively. 
The band-filling is a quarter in terms of holes. The pa- 
rameter to denotes the bare transfer integral, which is 
assumed to be uniform for simplicity, U the on-site re- 
pulsion strength, a and /3 the e-ph coupling strengths of 
modulating the transfer integrals by the molecular dis- 
placements (bond-coupled phonons) and the site energies 
by the molecular deformations (charge-coupled phonons), 
respectively, Ka and fT/j the corresponding spring con- 
stants, and uja and lu/s the corresponding bare phonon 
frequencies. In reality, the transfer integrals would be 
modulated by the molecular deformations as well, but 
such details are ignored. The static properties of this 
model {oJa=(jJi3=0) have already been studied so inten- 
sively, including exact-diagonalizationi^iM and density- 
matrix-renormalization-group^ studies. For charge or- 
ders in quarter-filled organic compounds, see for example 
Ref. [H. 

In the (EDO-TTF)2PF6 salt, displacements of an- 
ions in the insulator phase are substantial^^ and favor 
the (0110) (called bond-charge-density wave in Ref. [Tsl) 
ground state because the neighboring site energies are 
modulated in phase. This effect of anions counterbal- 
ances that of the nearest-neighbor repulsion favoring the 
(1010) (called Akp charge-density wave spin-Peierls in 
Ref. [ij) ground state. We here ignore such details for 
simplicity and assume that the effect of anions would 
renormalize the parameters related to the charge-coupled 
phonons, f3, K^, and top. Then, most of the model 
parameters can be deduced from the charge dispropor- 
tionation and the optical conductivity spectrum, as per- 
formed in Ref. 10: to=0.16, C/=0.93, a^/Ka=O.U, and 
P'^ /K/3=0.55, in units of eV. The classically treated, 
bond- and charge-coupled phonons are obtained by im- 
posing the Hellmann-Fcynman theorem. The opti- 
cal conductivity spectrum is then obtained with peak- 
broadening parameter set at 0.1 and shown in Fig. [Tfa). 
Because the on-site repulsion U is assumed to be the 
largest energy, the high-energy peak corresponds to the 
charge-transfer excitation to the (0200) configuration 
(wo20o ^^1-4), while the low-energy peak to that to the 
(1010) configuration (wioio ^0.4). 

III. ADIABATIC POTENTIALS 

The spring constants and i^T/j are so chosen that 
Uj I and I Vj \ are of the order of unity: -STq =0.035 
and Ki3^0.1 to give Uj=0.59, 0.85, -0.85, -0.59, • • ■ 
and Wj=-1.01, 1.01, 1.01, -1.01, • ■ • for the 12-site pe- 
riodic chain. The charge density per site is 0.07, 0.93, 
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FIG. 1; (Color online) (a) Optical conductivity in (0110) 
ground state with 12 sites. Schematic charge distribution (cir- 
cles) and deformation-induced potential (line), for (b) (0110) 
ground state, (c) (1010) photoexcited state, and (d) (0200) 
photoexcited state, (e), (f), and (g) are their adiabatic po- 
tentials as a function of amplitudes of bond-(uo) and charge- 
(«o) coupled phonons. 



0.93, 0.07, • • • in the ground state, though it is simply 
denoted by (0110). The charge distribution (rij) and 
the deformation-induced potential —f3vj are schemati- 
cally shown by the circles and the line, respectively, in 
Figs. [Ub)-md) for the ground and photoexcited states. 
The corresponding adiabatic potentials are plotted in 
Figs.[T](e)-mg). They are calculated with varying Uj and 
Vj through similar extension/contraction from the opti- 
mized configurations mentioned above. The maximum 
values of Uj and Vj are denoted by uq and vq, respec- 
tively. 

The minimum in the adiabatic potential of the (1010) 
photoexcited state is located at the origin [Fig.[T](f)] . The 
excitation energy is lowered by delocalizing the photocar- 
riers, so that both of Uj and Vj are strongly suppressed. 
In contrast, the minimum in the adiabatic potential of 
the (0200) photoexcited state is shifted from that of the 
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ground state to reduce uj only [Fig. [TJg)]. The energy 
contribution from the deformation-induced potential is 
basically unchanged by the charge-transfer process shown 
in Fig. [IJd). Meanwhile, the displacements shortening 
the distance between the neighboring hole-rich sites and 
stabilizing the tetramers are relaxed. Therefore, the shift 
of the stable phonon configuration is quite anisotropic by 
the (0200) photoexcitation. 



IV. PHOTOINDUCED MELTING DYNAMICS 
A. Photoexcitation by an oscillating field 

Photoexcitations are introduced by adding the mod- 
ified transfer term with the Peierls phase and the cou- 
pling with the molecular displacements (bond-coupled 
phonons) to the Hamiltonian, 



-^[to-a {uj+i - Uj)] iexp ^ / dtE{t) 
xc]_^Cj+i,o- + h.c. - ^ e{nj)5ujE(t) , 
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(2) 



where e is the absolute value of the electronic charge, a is 
the lattice spacing, 5=0. 05a is introduced to recover the 
length dimension. The time-dependent electric field E{t) 
is given by E{t) — — _Eoxt sinwoxt^ with amplitude E'oxt 
and frequency Wext for < t < Ti„ [E{t) is zero other- 
wise], where the pulse width is set at Tin- = 27riVoxt /woxt 
with integer iVoxt- 

It is easily shown by the gauge transformation that 
the vector potential introduced above is equivalent to the 
scalar potential.— The vector potential is convenient be- 
cause the boundary condition remains periodic. Of the 
two terms in Eq. ^ , the first is much larger than the sec- 
ond. Even if the second term is omitted, the numerical re- 
sults below are almost unchanged. The pulse width Tji-r is 
set to nearly coincide with the experimentally used value 
of about 200 in the present unit. The time-dependent 
Schrodinger equation for the exact many-electron wave 
function is numerically solved by expanding the expo- 
nential evolution operator with time slice dt=0.01 to the 
15th order and by checking the conservation of the norm 
and of the total energy for t > Ti„. The classical equa- 
tion of phonon motion is solved by the leapfrog method, 
where the force is derived from the Hellmann-Feynman 
theorem. The phonon frequencies lUq, and io^ are varied 
in such a way that one of the renormalized frequencies is 
close to the experimentally observed value of about 0.01. 
(either Wq or is set to be 0.02.) 



B. Numerical results 

First, the (0110) ground state is excited with a'oxt=l-5 
and A*'ext=50. The evolution of charge densities at 
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FIG. 2; (Color online) Time dependence of charge densi- 
ties, (a), (b), and (c) on short, intermediate, and long time 
scales, respectively, for excitation energy ajcxt = 1.5 close to 
'i'0200, A'oxt=50, and _Eext=0.340 just below threshold level, 
for u!ci=(^0~O.O2. 



an early stage is shown in Fig. [H^a). The period of 
E{t) is 27r/a;ext —4.2, which is close to that of the 
(0200) excitation, 27r/ti;o2oo —4.6, and far from that of 
the (1010) excitation, 27r/ti;ioio ~15. The bare trans- 
fer integral has the time scale of 27r/io —39. Because 
of Wext ^0200 J the charge densities at the hole-rich 
sites are largely oscillating with period 27r/wcxt, like 
(0110)^(0200)^(0110)^(0020). The oscillation at the 
hole-poor sites is initially small. This forced oscillation 
mainly at the hole-rich sites continues until the field is 
turned off at t=Ti„ ~210 [Fig. [2]^b)]. By this time, the 
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charge-density difference between tfie fiole-ricli and the 
hole-poor sites is reduced by the weakened Uj and Vj. 
Even for t > Tiir, the rapid oscillations due to charge- 
transfer excitations are evident. In addition, the slow 
oscillation on the phonon time scale, Tphonon —600, is 
clearly seen in Fig. [2ljc) . This coherent oscillation is al- 
lowed by the strong coupling between the charge and 
phonon dynamics. Here, -Ecxt is chosen to be just below 
the threshold for melting. Then, after several periods of 
the phonon oscillation, the charge configuration is nearly 
disordered and the charge order is quite reduced. 

How much the charge order is reduced can be quan- 
tified by averaging the charge density at the (initially) 
hole-poor sites over the period of the phonon oscillation 
at t <4000 (corresponding to an interval inside 2 to 3 ps) 
and by subtracting the initial density from it. It is plotted 
in Fig.[3]as a function of the increment in the total energy 
AE, for different excitation energies Woxt [indicated by 
the arrows in Fig.[l](a)] and different phonon frequencies 
LUa and LUf}. The arrows here show the energy difference 
between the undistorted and the ground states, which 
gives the lower limit of energy required for the melting. 
In all cases, the responses are linear for small AE, and 
their linear coefficients are nearly the same. For substan- 
tial AE, the curves are convex, indicating cooperativity 
in the photoinduced dynamics. It should be noted that, 
even if the charge order is almost completely destroyed by 
averaging over the phonon oscillation period, snapshots 
of charge densities and amplitudes of bond- and charge- 
coupled phonons show quite disproportionate configura- 
tions. Thus, the photoinduced melting does not produce 
a regular metallic phase. 

For excitations of energy cJcxt^O.S close to ujiqiq, the 
curve is insensitive to the combination of phonon fre- 
quencies [Fig. [2)Ja)]. Compared with the other excitation 
energies discussed below, the charge order is melted by 
smaller AE (0.31 < AE <0.35). In other words, the ef- 
ficiency of photoinduced melting is the highest on aver- 
age. This fact is similar to the fact, in one-dimensional 
half-filled dimerized Mott insulators, interdimer charge- 
transfer excitations reduce the spin-Peierls order much 
more efficiently than intradimer ones, leading to a pho- 
toinduced inverse spin-Peierls transition.^* In the present 
quarter-filled case also, carriers photogenerated by (1010) 
charge-transfer excitations are delocalized, so that the 
adiabatic potential has a minimum at the undistorted 
phonon configuration [Fig. [ijf)] to maximize the kinetic 
energy gain. 

For excitations of energy Ci;oxt=0.8 on the high-energy 
side of the (1010) charge-transfer peak, the efficiency of 
photoinduced melting is lowered and depends weakly on 
the phonon frequencies [Fig. [3I^b)]. The lower efficiency 
than the case with Woxt closer to wioio is reasonable in the 
sense that off-resonant excitations contain processes not 
simply represented by the (0110)^(1010) charge trans- 
fer: extra energy is needed to transfer charge. The de- 
pendence of the efficiency on the phonon frequencies has 
some regularity. Here, the curves for uja > are plotted 
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FIG. 3: (Color online) Reduction in charge order averaged 
over the period of oscillation at t <4000, as a function of incre- 
ment in total energy, for (a) Wext=0.5, A'ext = 17, (b) ^i;oxt=0.8, 
Ncxt='27, and (c) Li;cxt = 1.5, A^cxt=50. Bare phonon frequen- 
cies are (u^, tj;3) = (0.02, 0.02) (black), (0.02, 0.005) (red), 
(0.06, 0.02) (magenta), (0.005, 0.02) (green), and (0.02, 0.06) 
(blue) . The arrows show the lower limit required for the melt- 
ing explained in the text. 



in warm colors, while those for cua < i^fs in cold colors. 
The former case generally has a higher efficiency than the 
latter. 

This fact is more conspicuous for excitations of energy 
Wext=l-5 close to a;o20o [Fig-El^)]- Among the three exci- 
tation energies, the charge order is melted by the largest 
AE on average, indicating the lowest efficiency. However, 
the curve strongly depends on the phonon frequencies, so 
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that the average efficiency is not so meaningful. In fact, 
for (wq,, ujp)={0.02, 0.005), the efficiency is higher than 
the corresponding case of LUcxt—0-8- Meanwhile, the effi- 
ciency for {uja, 'j-'/3) = (0.02, 0.06) is very low. The origin of 
such dependence on the phonon frequencies is discussed 
below on the basis of the adiabatic potential, although 
the phonons do not necessarily evolve in its steepest de- 
scent direction. 

As mentioned before, the adiabatic potential of the 
(0200) photoexcited state [Fig. [Hg)] has a minimum 
shifted from that of the ground state [Fig.flje)] to reduce 
Uj only. Thus, the initial force reduces Uj. For Wq, > ujp, 
the energy supplied from the field is first transferred to 
the faster Uj motion, and then to the slower Vj motion, 
reducing both phonon amplitudes and the charge order 
in an efficient manner. Otherwise, the initially induced 
Uj motion is slow, and the secondary Vj motion requires 
higher energy than the first one, which is inefficient. 

The experiment on (EDO-TTF)2PF6 shows a corre- 
spondence of the frequency of the coherent oscillation 
with that of a Raman-active optical mode,- but its mode 
assignment is not done yet. This is why we vary the 
phonon frequencies in the calculations. At first glance, 
it may seem reasonable to assume that the displacement 
(uj ) is generally more collective and has a lower frequency 
Ua than the deformation (vj ) regarded as an intramolec- 
ular mode. However, the situation is not so simple be- 
cause the displacements of anions contribute to stabilize 
the (0110) charge order, as mentioned earlier. 

In order to study the effects of anion displacements and 
nearest-neighbor repulsion on the phase diagram, we add 
to the model. 
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(3) 

where wi is the displacement of the Ith anion from equi- 
librium, 7 the corresponding e-ph coupling strength, 
the corresponding spring constant, and V the nearest- 
neighbor repulsion strength. We use the relations U — 
F=0.93 and (3'^/Kf3 + 27^/7^^=0.55 so as to minimize 
the modification of the optical conductivity spectrum, al- 
though V still sharpens and heightens the lowest-energy 
peak. The phase diagram is divided into the present 
(0110) order approximately for V < 0.1 + 7/1.3 and the 
(1010) order for V > 0.1 + 7/I.3. The almost linear 
relation is due to the competition between 7 favoring in- 
phase modulation of {1121-1) and (ri2i) and V favoring 
out-of-phase modulation of them. The strong reduction 
of the critical V is due to the strong Holstein coupling 
13'^ /Kp. Without e-ph couplings (a=/3=0), the (1010) 
order needs V > 2tQ in the limit of large U and a some- 
what larger V for finite U. Therefore, without anion 
displacements (7=0), the nearest-neighbor repulsion of 
only 10% of the magnitude of U is thus found to convert 
the (0110) order into the (1010) order. Because the an- 
ion displacements strengthen the e-ph coupling /3, they 
make the Vj phonon effectively heavier. Thus, it is pos- 



FIG. 4: (Color online) Time dependence of charge densities 
for uja~0.02 and uj0—O.O6. Photoexcitations are (a) ajcxt=0.5 
close to tJioio, A'^ext^lT, and £0x1=0. 030; and (b) a;oxt=1.5 
close to CJ0200, Ai'cxt=50, and iJoxt=0.280. 



sible that the Vj phonon has a lower frequency uj/s. If 
it is the case, the efficiency with LUcxt—1-5 can be higher 
than that with ijjoxt=0.8. It should be noted that, if the 
reduction in the charge order is plotted as a function of 
the number of absorbed photons (AiJ/wcxt), the excita- 
tions with Wcxt = l-5 are the most efficient and those with 
are the least efficient on average. 

Details of the actual dynamics of charge densities and 
phonon amplitudes are more complicated than those 
naively expected from the adiabatic potential as de- 
scribed above. The charge dynamics for Wq, ~ LOfs were 
close to a sinusoidal evolution [Fig. [H^c)], but they are 
not always the case. Below we will describe the actual 
dynamics for to a ^fi- 

For uja < 1^/3, the charge dynamics sensitively depend 
on the excitation energy (Fig. |4]). When the excitation 
energy is close to wioioi 'j-'cxt=0.5, the main component 
in the charge-density oscillation is due to the Vj phonon 
with the higher energy [Fig. H^a)]. Because the corre- 
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FIG. 5: (Color online) Time dependence of charge densities 
for uja=0.06 and up— 0.02. Photoexcitations are (a) ajcxt=0.5 
close to ojioio, A''ext = lT, and i5cxt=0.032; and (b) ajcxt=l-5 
close to CJ0200, A''ext=50, and iJoxt=0.310. 



sponding adiabatic potential is rather isotropic with re- 
spect to the Uj and Vj ampUtudes [Fig.[Hf)], the phonon 
with a stronger couphng strength > a'^/Ka), 

Vj, more contributes to the charge-density oscillation. 
On the other hand, when the excitation energy is close 
to wo200j '^cxt=l-5, the main component in the charge- 
density oscillation is due to the Uj phonon with the lower 
energy [Fig. HJb)]. This is because the minimum in the 
corresponding adiabatic potential is shifted from that of 
the ground state to reduce uj only [Fig. W[s)] ■ 

For LOa > i^f3, the charge dynamics are much more 
insensitive to the excitation energy (Fig. [5]) than for 
Ua < described above. The charge-density oscil- 
lation is dominated by a single component, whose fre- 
quency weakly depends on tJoxt- The frequency of this 
coherent oscillation for tJoxt=0.5 [Fig. [5l^a)] is slightly 
lower than that for tJoxt^l-S [Fig. (Hlb)]. This small dif- 
ference would be caused by the difference in the cou- 
pling strengths, 0^ /Kp > a^/Ka, and the difference in 



the shapes of the adiabatic potentials, as in the case of 
LUa < (^0. The charge dynamics for (cJq,, a;;3)=(0.02, 
0.005) are slower but otherwise similar to the present 
case for {uja, uj(3)~{0.06, 0.02). Such (wq,, a;^)-dependent 
theoretical dynamics would be helpful in assigning exper- 
imentally observed oscillations. 



V. SUMMARY 

Photoinduced melting of charge order is calculated ex- 
actly for electrons and classically for phonons in the one- 
dimensional quarter-filled Hubbard model with different 
types of e-ph couplings. The strengths of the e-e and 
e-ph interactions and the transfer integral are taken to 
reproduce the charge disproportionation and the opti- 
cal conductivity in the (0110) ground state of (EDO- 
TTF)2PF6. The efficiency of (0200) photoexcitations is 
generally (but not always) lower than that of (1010) pho- 
toexcitations, which generate delocalized carriers. Rel- 
ative dynamics of the two phonon amplitudes are influ- 
enced by the shape of the adiabatic potential and the rel- 
ative phonon frequency. As a consequence, the efficiency 
of (0200) photoexcitations can be high if the frequency 
of the bond-coupled phonons is sufficiently higher than 
that of the charge- coupled phonons. On the other hand, 
the efficiency of (1010) photoexcitations is insensitive to 
the relative phonon frequency, though it is lowered if the 
excitation energy is off- resonant. 

Here, we show that the relative phonon frequency can 
become an important parameter characterizing the pho- 
toinduced melting dynamics if the order is stabilized by 
different types of phonons. The charge-density oscilla- 
tion is also affected by the relative phonon frequency, so 
that theoretical calculations are generally useful for the 
mode assignment of the coherent oscillation and at least 
for the identification of interactions stabilizing the order. 
Furthermore, consideration of relative phonon frequen- 
cies may open the possibility for controlling the efficiency 
of photoinduced dynamics. If many types of phonons 
are involved with photoinduced dynamics, we conjecture 
that efficient structural deformation is achieved when the 
phonon coupled with the initial charge-transfer process 
has the highest frequency, the secondary phonon induced 
after the first has the second highest frequency, and so 
on. 
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